Radiation from Excited Vortex in the Abelian 



Excitation of a vortex in the Abelian Higgs model is investigated with the help 
of a polynomial approximation. The excitation can be regarded as a longitudinal 
component of the vector field trapped by the vortex. The energy and profile of the 
excitation are found. Back-reaction of the excitation on the vortex is calculated in 
the small k limit. It turns out that in the presence of the excitation the vortex 
effectively becomes much wider - its radius oscillates in time and for all times it 
is not smaller than the radius of the unexcited vortex. Moreover, we find that the 
vector field of the excited vortex has long range radiative component. Bound on the 
amplitude of the excitation is also found. 
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1 Introduction 



There are many reasons for detailed studies of dynamics of vortices. Let us mention here 
the cosmic string hypothesis IJ, condensed matter physics [0, QCD flux-tube || and a 
search for a physical model of a relativistic string Q. 

Cumulative efforts by various groups of physicists have unravelled many aspects of 
fascinating dynamics of vortices. Vortex is a quite complex object - due to its macroscopic 
thread-like extension it has nontrivial time-evolution with occasional self-interactions of 
various kinds like scattering || or interconnection ||. Moreover, because of its finite 
transverse size it can also act as effective potential well trapping various fields. The first 
example of this kind was the superconducting string . It turns out that this can happen 
also in the more restrictive framework of Abelian Higgs model ||, |9|. Here the trapped 
fields are just some components of the vector field. Because other components of the same 
vector field form the vortex itself, it is natural to regard this phenomenon as the excitation 
of the vortex. This aspect of dynamics of the vortex becomes important when the vortex 
is subjected to such violent transitions like the interconnection or hard scattering on an 
external potential. In fact, numerical investigations of the interconnection clearly show 
that after this process the involved parts of vortices are excited, [ lOfl . The vortex can also 
become excited by interaction with particles scattered on it. 

Excitations of the vortex in the Abelian Higgs model were studied recently in the 
paper [11]. In that paper an effective membrane description of the vortex was proposed. 
The membrane was identified with the surface on which the energy density is maximal - 
for a straight-linear, static vortex with winding number \n\ > 2 it is a cylinder of finite 
radius. In the membrane model there exists solution which describes a symmetric bulb 
travelling along the straight-linear cylinder. It corresponds to a local excitation of the 
vortex. The shortcoming of such effective model is that it does not yield information on 
detailed form of the excited vortex fields. 

Comprehensive numerical study of various excitations of local and global vortices has 



been carried out by the authors of Ref . [ 12 ] . The excitation we consider in the present 
paper has not been investigated there. 

In general, excitations of the vortex can be divided into two classes, according to the 
involved components of fields. In the first class there are excitations which do not involve 
new components of the gauge or Higgs fields - only modes of those components which 
form the vortex itself are excited. For the vortex these are the modulus of the Higgs field 
and the Ag component of the vector field (8 denotes the azimuthal angle). The other class 
of excitations requires to extend the vortex Ansatz by allowing for non-trivial values of 
additional components of the Higgs or the vector fields, i.e. the phase of the Higgs field 
or the A z component of the vector field (the A component is equivalent to the phase 
by a gauge transformation). These latter excitations should rather be regarded as bound 
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states of the additional fields with the vortex. The excitations are usually investigated in 
a linear approximation - they are treated as small corrections to the basic vortex fields. 

The excitation we are interested in belongs to the second class - it involves the A z 
component. It can be regarded as a bound state of longitudinal component of the vector 
field with the vortex. In the paper || the existence of the excitation has been proved, 
and it has been shown that it can be represented by a 2-dimensional Proca field which 
can propagate along the vortex. These results have been obtained in a linear approxi- 
mation with respect to deviations of the fields from the static vortex solution. Merely a 
rough description of the excitation was given - the main obstacle was the lack of analytic 
expression for the vortex solution even in the Bogomol'nyi limit [13]. Many interesting 
questions were left without answer, e.g., what are properties of the excitation beyond the 
linear approximation. 

The excitation seems to exist for all values of the parameters of the Abelian Higgs 
model compatible with the existence of the stable vortex of topological charge +1 |§. 
Nevertheless, in the present paper we concentrate mainly on the case of very small k, where 
it is possible to give explicit and relatively simple analytic formulae. Here k = mA_/mH, 
where uia and m# denote the masses of the vector and Higgs particles, respectively. 

We would like to present a detailed study of the excitation in the case k < 1/2. Our 
considerations are based on approximate, analytic expression for the vortex functions 
obtained in (appropriately adapted) so called polynomial approximation, used till now for 
domain walls |14j]. Next, we obtain formulae for the profile and energy of the excitation. 
This in turn makes it possible to calculate back-reaction of the excitation on the vortex. 
The back-reaction is calculated in the case of very small k when certain simplifications 
occur. We find that in the presence of the excitation the vortex has oscillating radius, 
and we calculate the frequency and amplitude of these oscillations. We also notice that 
there seems to exist an upper bound on amplitude of the excitation. Such a bound reflects 
the non-linear character of the vortex dynamics and it could not be found in the linear 
approximation. However, the most striking result is that the excited vortex field has 
a long range component which is interpreted as radiation of the vector field. Also this 
phenomenon is absent in the linear approximation. 

The general conclusion of our paper is that the effects of the back-reaction are very 
important - they qualitatively change the picture of the excited vortex obtained in the 
linear approximation. 

The plan of our paper is as follows. In the next Section we present the polynomial 
approximation for the vortex. In Section 3 the profile and energy of the excitation are 
calculated. In Section 4 we investigate the back-reaction of the excitation on the vortex. 
In Section 5 we have collected several remarks about the method used in our paper and 
about the obtained results. 
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2 The polynomial approximation for the unexcited 
vortex 



We consider Euler-Lagrange equations of the Abelian Higgs model: 

A 2m 2 

{8 V + iqA v W + iqA v )* + -$(|$ | 2 - — ) = 0, (1) 

= iq(<S>*d v <S> - $9"$*) - 2q 2 A"\<5>\ 2 , (2) 

where m,q,\ > 0. The mass of the Higgs particle is equal to m 2 H = 2m 2 , while for 
the vector field m\ = K 2 m 2 H) where k 2 = 2q 2 /\. Signature of the space-time metric is 
(+1,-1,-1,-1)- 

The solution of Eqs. (1),(2), representing an infinite, straight-linear, topological charge 
+ 1 vortex lying along the x 3 -axis, is sought with the help of the axially symmetric Ansatz 



* = ^e»F{p\ (3) 
A = A 3 = 0, (4) 

A^-H(p), A 2 = --H(p), (5) 
P P 

where p = yj (x 1 ) 2 + (x 2 ) 2 is the radius in the (x^x 2 ) plane and 6 = arctan^/x 1 ) is the 
azimuthal angle. The functions F, H obey the following boundary conditions: 

F(0) = if (0) = 0, (6) 

F(oo) = 1, lim pH(p) = --. (7) 
From (1-5) we obtain equations for F and for \ = QpH(p) + 1: 

F" + y-^X 2 F+ 1 -(F-F 3 ) = 0, (8) 

X"---K 2 F 2 X = 0, (9) 
r 

where r = m H p, and ' denotes d/dr. Boundary conditions for \ have the form 

X (0) = 1, x(oo) = 0. (10) 
The magnetic field of the vortex is B l = 5 l3 B(r), where 

B = (11) 
q r 
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The main idea of the polynomial approximation is to approximate the fields inside 
the axially symmetric vortex by simple polynomials in r, and to match the polynomials 
smoothly with asymptotic form of the vortex fields found for r > tq, where tq is the 
matching point. We shall call r the matching radius of the vortex. Accuracy of such 
approximation in general increases with the order of the polynomials. In the present 
Section we shall assume a third order polynomial for F and a sixth order polynomial for 
X- Inserting these polynomials in Eqs. (8), (9) we find that the polynomial for F contains 
only odd powers of r, and the polynomial for x only even ones. Thus, denoting these 
polynomials by underlined letters, we have 

E = /ir-|/ 3 r 3 , (12) 

X = l-lh 2 r 2 + ^h 4 r 4 -^h 6 r 6 . (13) 

The order of the polynomial x is th e maximal one compatible with the assumed third 
order for the polynomial F_ — in order to compute the term ~ r 8 in x we would have to 
know the function F in Eq. (9) up to r 5 . Equations (8), (9) give the following recurrence 
relations for the coefficients of the polynomials 

h = \{\ + h 2 )h, 

K = 3k 2 / 2 , (14) 
he = 5« 2 / 1 (2/ 3 + 3/ i2 / 1 ) 

The asymptotic values F, x of the fields F, x outside of the vortex, i.e. for r > ro, where 
r will be determined later on, we find using Eqs. (8), (9) and the boundary conditions 
(7), (10). There is a complication due to the fact that the asymptotics of the Higgs field 
depends on k ||17|| . For k < 1/2 it is governed by the x 2 F/r 2 term in Eq.(8), while for 
k > 1/2 it is governed by the Higgs field mass term. In the following we shall assume that 
k < 1/2. This is mainly to limit the body of our paper - our methods of calculations of the 
properties of the excitation can be applied for any k. An additional reason is that for small 
k certain simplifications are possible, e.g., one can use a kind of adiabatic approximation. 

If for the Higgs field we take 

F (0) = 1, (15) 
then in the asymptotic region Eq.(9) is effectively simplified to 

x ^x = o, 

r 

and its vanishing for large r solution is given by 

X (0) = CorK^nr), (16) 
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where K\ is the modified Hankel function [15] and Cq is a constant (to be determined 



later on). In the next step, one could insert the asymptotics (16) for x m t° Eq.(8) and 
to determine corrections to the leading asymptotics (15). Iterating these steps and taking 
higher order polynomials for F_ and x we can obtain a very accurate description of the 



vortex which will be presented in a separate paper |L6f . In the present paper, devoted to 
the excitation rather than to the vortex, we will be satisfied with rather crude description 
of the unexcited vortex. 

The last step in obtaining the static vortex solution is to match the asymptotics with 
the polynomials (12), (13). We know from []I7| that for the exact static vortex solution F is 



real analytic in (x 1 , x 2 ). As for the magnetic field, it is also true at least in the Bogomol'nyi 
limit. For this reason, the number of the matching relations depends only on the number 
of arbitrary constants present in the asymptotic solutions and in the polynomials - we just 
require maximal possible smoothness of the fields within the adopted Ansatz (12,13,15,16). 
In the present paper we will approximate F and x by functions of r of the C 1 class at 
r = ro, and real analytic for r ^ ro- More accurate and smoother approximations to the 
vortex fields will be presented in [|16| . Thus, we choose the following matching conditions 
for F and x : 

F(ro-) = F(r +), F'(r -) = F'(r„+), (17) 

X(ro-) = X(r +), X(r ~) = x'{r +). (18) 

Here g(r ±) = lim r ^ ro ± g(r) with g = F,x- The radius r and the coefficients Co, fi, hi 
are to be determined from (17), (18). Conditions (17) with (15) taken into acount give 

A=2P /3 = ^- (19) 

F<® = *L-±(L)\ (20) 
2r 2W 1 ^ 

Next, we calculate h^h^ and h§ from formulae (14) - in this way we find 



I I 4 \ 2 9 h 2 4 3 k 2 (I 1 . , 



The final step is to obey the matching conditions (18) for x- Taking x = X^ we fi n d 
that they are equivalent to the following two equations from which we can determine cq 
and ro 

~ I + I (1 + i^ 2)r ° + lls^ = c ° r o^i( Kr o), (22) 



8 - ¥ 2 - 2 £n 2 r* K>(kt 
l-!(l + |^)r 2 -iis« 2 r 4 ^ + K 1 {kt 



1 + Wq—- r. (23) 
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Equation (23) can be easily solved for tq by numerical methods for any given k. For very 
small k one can also give analytic formula: 

1 fi 

rl = ^-{l + K 2 \nK) + 0{K 2 ). (24) 

(we have used the formula K±(z) = l/z for < z « 1). Thus, for k — > the Higgs radius 
of the vortex reaches the fixed value 4/v / 3 ~ 2.31. Comparison with numerical solutions 
of Eq.(23) shows that formula (24) gives slightly too big value of tq, and that the error is 
only 6% even for k = 0.3. For k — 0.1 the error is 1% and it decreases rapidly for smaller 
k. In general, the radius decreases with increasing k, see Fig.l. 




K 



0.1 0.2 0.3 0.4 0.5 



Fig.l. r as a function of k, determined from Eq.(23). 

The constant c is found from Eq(22). In particular, for k — > 

4 

c = k(1 - -^k 2 \yik) + 0{k 3 ). (25) 
v3 

The magnetic field is given on outside of the vortex by formula (11) with \ — 
(formula (16). It has indefinitely increasing range in the limit k — > 0. At the same time 
its magnitude decreases so that its integral gives the unit of the magnetic flux. 

Let us remark that for k < 1/2 a more accurate approximation for the asymptotic 
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field F is provided by the formula 
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Using this asymptotics instead of F^ in the matching conditions (17), (18) and a fifth 
order polynomial for F_ we find that for k — > r = 2.11. Thus, the improved asymptotics 
gives slightly smaller matching radius. Such more accurate asymptotics F and a fifth 
order polynomial F_ are used in Section 4. Yet more accurate asymptotics is given by the 
formula 



2 



£(0)2 + cKq{ 



r) 



where c is a constant. This asymptotics is considered in [16 | . 

The approximate, analytic form of the vortex functions F, x we have obtained in this 
Section for k < 1/2 will be used to study the excitation in more detail than it was done 
in m. 



3 Analytic description of the excitation in linear ap- 
proximation 

The excitation is another vortex-like solution of Eqs. (1,2). It belongs to the same topo- 
logical class as the vortex itself. Arguments for existence of such a solution, based on a 
linear approximation, are given in |§,|J. This solution is time-dependent. 

The Ansatz (3-5) is now extended by including A 3 = A(t, r), while Aq = again. The 
F and % functions are now time-dependent. In this case equations (1),(2) reduce to the 
following set of equations 

^F + F" + — - ( ^ X 2 + -V A 2 \F+ ^-(F-F 3 ) = 0, (26) 



r 




-k 2 F\ = 0, (27) 
- k 2 F 2 A = 0, (28) 

where t = mnx and the dots denote d/dt. Notice that there is no direct coupling of the 
A field to the x field, and that A influences F only by the single term ■^ r A 2 F in Eq.(26). 

We shall seek for the solutions of this set of equations using an iterative procedure. 
First, in the present Section we aproximately solve Eq.(28) assuming that the fields F and 
X form a background not influenced by the A field, i.e. that they are equal to the initial 
vortex fields discussed in the previous Section. With the back-reaction of the A field on 
the vortex switched off, equation (28) becomes linear equation for A with explicitely given 
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function F(r) - hence the name "linear approximation". The back-reaction of the A field 
on the F and \ fields we shall calculate in the next Section. We hope that if we continued 
the iterations we would obtain a sequence of approximations convergent towards the exact 
solution at least for sufficiently small -^p-A 2 (which is a dimensionless quantity). 

Thus, in the present Section we will find the bound-state type solution of the equation 

_ a + A" + — - k 2 F^ 2 A = 0, (29) 
r 

where F^ denotes the Higgs field of the unexcited vortex, calculated in the previous 
Section. We shall use the crude asymptotics (15) and the third order polynomial (12) for 
the Higgs field - it is sufficient in order to find the excitation and its basic characteristics. 
More accurate form of the Higgs field would require significantly longer calculations and 
we think that the improvement would yield merely more accurate numbers, with no change 
in the overall picture of the excited vortex. If r < r , the functions F, x and A are again 
approximated by polynomials F_, x and A, respectively. For F_ and x we have formulae 
(12,13), and 

A = a (t) - \a 2 {t)r 2 + ^a 4 (t)r 4 - ^a 6 (t)r 6 . (30) 

We have noted that the coefficients in A can be time-dependent. Inserting formula (30) 
in Eq. (28) we obtain the following set of equations 

d = -2a 2 , (31) 
4 

a 2 = --a 4 + 2K 2 a f 2 , (32) 
6 

a 4 = --a 6 + 8/t 2 a /i/ 3 + 12k 2 a 2 f 2 ■ (33) 

If r>r , in Eq.(29) we put F (0) = 1. It follows that A — A, where 

A — ci cos(u t + S)K (k r). (34) 

Here c\ is a constant, Kq is the modified Hankel function and 

ul + k 2 G = k 2 . (35) 

The frequency u and the constant c are determined from matching conditions at 
r = r : 

A(r -) = A(r +), (36) 
A'(ro-) = A'(r +). (37) 

We assume that the coefficients a 2 fc, k = 0, 1, 2, 3, have the same time-dependence as A, 
i.e. 

a 2 k = «2fc cos(a;ot + 5) 
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with constant a 2 k- Then, using Eqs. (31-33), (19) we can express a 2 , a^a Q by a ,r Q ,K. 
We obtain that 

A = a(r) cos(co t + 5), (38) 

where 



a(r) = a 



1 2 2 1 

4^ + 64 
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2 I 4 

k r 



'45 UqK 2 



48 I 8 r, 



12- 



1 



6 I Jo 







+ g^n I r 



The matching conditions (36), (37) are equivalent to the following two equations 

Cl 



1 1 2 2 , 1 4 4 ^ 

1 - 4^° + 64^° - 384 



fi fi ^22 224 

UJ ° r ° ~Q4 Kr °~ 128^°* r ° 



a 



(39) 



K (k r ), (40) 



2^2 



1, ,4„4 



45, ,2^2„4 



0' 



15 ,^2^2 



8^0' ~ 32^0' 



32^0' 







nr. 



o 



1,,,2 2 I J_,,,4 r 4 _ 

A uj r -+- m w r 



1 , S r 6 _ 7 2^2 _ 15 ,2^2^,4 
SSi^O'O 64 128^0^ '0 



-2r k 



K' (k r ) 
^o(Vo)' 



(41) 



where K' Q (z) = ^K (z) with z = A; r , and r is to be calculated from Eq.(23) of the 
previous Section. 

From Eq.(40) we can determine only the ratio c\/ao, so the amplitude «o remains 
arbitrary in accordance with linearity of Eq.(29) with the fixed F^°\ 

Equation (41) gives ooq as a function of k. In general, one has to use numerical methods 
to find this function. It is plotted in Fig.2. 




Fig.2. oj as the function of k, determined from Eqs. (41), (35) 
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If k <C 1, one can use formula (24) and the following asymptotic formula for the 
modified Hankel function |lT5fl for < z << 1 



zK' (z) 



K (z) ln( 7 z)-ln2' 

where 7 ~ 1.781072 is the Euler constant. Formula (35) implies that for k — > ujq and ko 
vanish. Keeping only the leading terms on both sides of Eq.(41) we obtain the equation 



: k + -k )r 



v 16 2 u/ u logfc ' 
which gives 

~ exp 



32 1 1 1 

'23^2^/ 

It is easy to see that the value of the coefficient in front of this formula for ko depends on 
subleading terms in Eq.(41). Taking these terms into account we find that for k<C 1 

,9 1 f 32 1 1 1 

kl = 3.86^ exp ---- , 42 

?0 [ 23 K Tq J 

where r is given by formula (24). Comparison with numerical solutions of Eq.(41) shows 
that formula (42) gives slightly too big value of k , with 13% accuracy even for k = 0.3. 
For k — 0.1 the error is 2% and it decreases rapidly for smaller k. 
From Eq.(40) we find that for small k 

- = -f^ 2 - («) 

a lb 

It follows that when k — > 0, the function ciK (ko,r) for any finite fixed r tends to a 
constant value equal to «o, hence A — > ao- 

The existence of the excitation in the linear approximation can be seen by regarding 
the equation 

1 , s 2 2 

1 / „ a\ K „/n^o OJn 



obtained from Eq. (29) by substituting A = cos(uj t + S)a(r), as Schrodinger equation 
for a particle in the 2-dimensional potential well given by V(r) = f F(°) 2 . It is a well 
known fact jl8| that any potential well on a plane has at least one bound state, hence 
localized non-zero solution A of Eq.(29). For small k the potential well has finite width 
r ~ 4/ a/3 and it is very shallow (V m3iK = k 2 /2). Nevertheless, the bound state exists 
and the corresponding frequency uq, is given by formula (42). Estimate made in 0] shows 
that the potential well is too weak to support p-wave bound states. Therefore, we do not 
consider more general, #-angle dependent A field. 
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4 The vortex in presence of the excitation 



Mathematically, the presence of the excitation changes the vortex because of the A 2 term 
in Eq.(26). Because A 2 ~ cos 2 (uj t + S) — |[1 + cos(2uj t + 25)}, the perturbation has two 
components: the static one and the one oscillating with the frequency 2u>o. In order to 
calculate their influence on the vortex we reanalyse Eqs.(26),(27) assuming that A is given 
by the approximate solution of Eq.(28) found in the previous Section. We shall consider 
the limiting; — > when we can do this analytically. 

In the case of very small k, A 2 changes in time very slowly as the frequency 2uj is 
very small. Because we expect that the induced by the presence of the excitation time 
dependence of F is also characterized by the frequency 2uj , the term F on the r.h.s 
of Eq.(26) is of the order lj 2 ~ k 2 and it can be neglected. Therefore, equation (26) 
approximately reduces to the following equation 

F" + -F' - (\ X 2 + ^rA 2 )F + \{F - F 3 ) = 0. (44) 
r r z m H 2 

Now let us turn to the Eq.(27) for the field x- Again, we expect that the induced time 
dependence of x is characterized by the frequency 2u)q, so that the term x is of the order 
k 2 . However, if we neglect this term then we also have to neglect the term k 2 F 2 x, and 
Eq.(27) reduces to the following equation 

X " - -X' = 0. (45) 
r 

Its general solution has the form 

X = c 2 + d 2 r 2 , (46) 

where c 2 , d 2 are arbitrary constants. It is easy to see that such x can n °t obey the boundary 
conditions (10). Thus, neglecting the terms ~ k 2 in Eq.(27) altogether is an oversimplifi- 
cation if we consider the whole range of r, i.e. [0, oo). On the other hand, it is acceptable 
approximation if we consider x in a finite interval of the r variable, e.g., in a vicinity of 
the matching point, now denoted by r e . It breaks down only if we try to apply it also for 
r — > oo. This means that the limits k — > and r — > oo are not interchangeable, and the 
correct order of them is that first we should find x f° r ver Y large r and only after that 
we can consider the limit k — > 0. One could see this by analysing the vortex solution of 
Section 2, cf. formulae (16), (25). Later on in this Section we shall find the asymptotics x 
for the very large r. For the current calculations it is sufficient to know that for r ~ r e x 
is given by formula (46). 

From Eq. (44) we obtain approximate form of the Higgs field F in the presence of the 
excitation. As always in the polynomial approximation we separately consider the regions 
of small and large r. The matching takes place at r = r e . In the region r < r e we will 
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approximate F by a fifth order polynomial, 

F = hr ~ ^/ 3 r 3 + ^hr\ (47) 

The reason is that the more precise (from that considered in Section 2) form of the asymp- 
totics of the Higgs field which is given below requires at least the fifth order polynomial if 
the matching conditions (17) are to be satisfied. Instead of the first of recurrence relations 
(14) we now have 



f *=l(l +h >-&/^^ < 48 > 

We also obtain the relation 



A = l { \ + fc2 " k al)h + ^ ,I4 + ^ + /? " ^ aott2)/l ' (49) 

Notice that passing to the fifth order polynomial does not increase the number of arbitrary 
constants in our Ansatz, because the new constant / 5 is related to the other constants by 
formula (49). 

For x inside the vortex we again use the polynomial (13), and now the coefficients h 2 k 
can in principle depend on time. However, the corresponding recurrence relations in the 
limit k — > reduce to 

h 4 0, he = 0. (50) 

The matching conditions (18) applied to the polynomial x an d X given by formula (46) 
imply that 

c 2 = 1, d 2 = ~^h 2 . (51) 
The asymptotics F of the Higgs field obeys the following equation 

F" + \P> - (lx 2 + ^)F + \{F - F 3 ) = 0, (52) 

where A is given by formula (34), and x is found below. We use the observation following 
from formulae (34), (43) that for any r > in the limit k — > 

(A 2 y ~ k 2 . 

Hence, A 2 changes very slowly with r for small k. As we shall see below, the same is true 
for Xi essentially because the mass of this field is equal to k. Then, it is easy to check that 
the equation for F has the following approximate solution 



F=Jl-- 2 x 2 -%A\ (53) 



in 



H 
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Adopting this solution is equivalent to neglecting the derivative terms (~ F', F") in 
Eq.(52). 

The next step is to obey the matching conditions (17). Because the matching takes 
place at finite r = r e we may use formula (46). We also use formulae (50), (51), and we 
notice that for finite r and k — > 

A 2 = a 2 , 

as follows from formulae (34), (42), (43). The matching conditions (17) give the following 
relation 

= 8(5^-12) 
,•2(32-^)^-2 



where 



as well as the equation 



0=1-^ + 2/12, 



m 



H 



5/3 - 6x ( 2 (3 (3 2 8 x 3 (5p - 6xf , 



for 



Equation (55) has the following solution 

x = do {3, 

where the constant do can be easily determined by numerical methods, 

d w 0.449. 

Moreover, we shall show later on that h 2 ~ 0. Hence, 

2.11 



It is clear from this formula and from the formula 

dp = al cos 2 (u t + 5) 

that r e oscillates in the interval [2.11, 2.11/(1 — 2q 2 a 2 l /m 2 H )]. Let us recall that 2.11 is the 
matching radius in the absence of the excitation, as mentioned at the end of Section 2. 

We see from formula (56) that as a 2 , approaches m 2 H /(2q 2 ) the amplitude of oscillations 
of r e increases indefinitely, and for still higher amplitude of the excitation formula (56) 



(56) 
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gives unphysical r e . We interpret this result as a hint that there is an upper bound 
~ fflj; / \yj2q) on the amplitude ao of the excitation. Such a bound would reflect nonlinear 
character of the set (26-28) - in the linear approximation no restriction on the amplitude 
appears. It would be desirable to check the existence of the bound by another method, 
e.g., by direct numerical analysis of Eqs. (26-28). Interesting question what happens to the 
vortex if the bound is exceeded we will leave for a separate investigation. 

In order to find the x field °f the excited vortex for the very large r we have to solve 
equation (27) with F = F, where F is given by formula (53). As discussed earlier, now 
we should not neglect the terms proportional to k 2 . Inserting formula (53) and neglecting 
the term 2/t 2 x 3 /r 2 (because we expect that x is small for r > r e in comparison with the 
other terms - subsequent computation confirms this expectation) we obtain the following 
equation 

-X + X"- -X' - m 2 effX + n 2 ^-K 2 (k r) cos(2^ t + 25)x = 0, (57) 
r m H 



where we have introduced the notation 



m 



q 2 c 2 



eff 



= K 2 (l-l^K 2 (k r)). (58) 



m 



H 



Next, we split x into a static and time-dependent parts, 

X = Xs + Xoj, (59) 
where the static part by definition obeys the equation 

Xs ~ \x s ~ m lffXs = 0. (60) 
It follows from formula (59) and Eqs. (57), (60) that Xui obeys the equation 
-L + XL - \Xu - m l ff Xu + K 2q -^K 2 (k r) cos(2^ t + 25)xu> 



= - K ^K 2 (k r) cos(2cu t + 25)xs- (61) 



Thus, Xs ac ts as a source for Xu- 

In the present paper we will discuss only a perturbative solution of Eq.(61), obtained 
by expanding Xui in powers of the dimensionless amplitude qao/rriH of the excitation. This 
solution is expected to be a good approximation to the exact solution of Eq.(61) if the 
dimensionless amplitude is small, so from now on we assume that this is the case. The 
perturbative solution is sufficient to show that the excited vortex radiates the x field. We 
start from the observation that the r.h.s. of Eq.(61) together with formula (43) suggest 
that Xui is of the order q 2 (x^jm 2 H . Therefore, the last term on the l.h.s. of Eq.(61) is of the 
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order (q 2 al/m 2 H ) 2 and can be neglected if the dimensionless amplitude of the excitation 
is small. For the same reason, m e ff in Eq.(61) can be replaced by k. Solving obtained 
in this way the simplified version of Eq.(61) we find the lowest order contribution to the 
perturbative series for Xw 

Let us introduce a new function h u (r) 

2 2 

X* = cos(2cu t + 26) r h u {r). (62) 



m 



The simplified form of Eq.(61) is equivalent to the following equation for the function h w : 

hi + -hi - \K + k\K = --^\K 2 (k r)x s (r), (63) 
r r z rag 

where 



k\ = \JAujq — k 2 . 

Observe that for the very small k the coefficient k\ is positive because then uj 2 ~ k 2 , see 
formulae (35), (42). 

The asymptotics x is determined by Xs and h u . Equation (60) has the following van- 
ishing at the infinity, approximate solution 

X s « csnrKi (m e ffr) , (64) 

where C3 is a constant to be determined later. 

General, vanishing at the infinity solution of Eq.(63) is given by the formula 

K{r) = c 3 /ioo(A;ir) + c A Ji(kir) + CsN x (kir), (65) 

where C4, C5 are constants to be determined later on; J%, Ni are Bessel and von Neumann 
functions, respectively f[5|] ; and h^z) (with z = k%r) is a particular solution of the 
inhomogeneous Bessel equation 

^# + + (i - -^(z) = -£4*i(£*)K2(r*)- (66) 

CtZ Z CLZ Z rC x Q^q tZ\ rC\ 

For small k we can approximate 



Equation (66) can be solved by a standard method [|19 



h 00 (z) = f(z)J 1 (z)+g(z)N 1 (z), (67) 

where 

\/3 k-\/3 
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M = ~t\ dxxN&W-^KSAx), (68) 
6«n ^2 v3 kV3 



and 

qrr^ /'CO q* U 

9{z) = ~d L dxxMx^i-^K^-^x). (69) 

The constants c 3 ,c 5 can be fixed by matching the asymptotics x with x, given by 
formulae (46), (51), at r = r e . This requirement gives the following equations 

Ure) + ^>r e h„(r e ) cos(2cu t + 25) = 1 - ^r e 2 , (70) 

X a (r e ) + ^-^[r e hi(r e ) + ^(r e )] cos(2^ t + 25) = -h 2 r e . (71) 
Next, we notice that for small k (when Kr e << 1) 

Xs(r e ) ~ c 3 , xl(r e ) ~ 0, 

2 

Kr e Ji(V3Kr e ) ~ 0, K,r e Ni(VSnr e ) t 



Then it is easy to see that conditions (70), (71) imply that 



/i 2 = 0, c 3 

The definition (69) gives for z = 



1-^, c 5 = -^(0)c 3 . (72) 



9(0) = (73) 



The constant C4 is not fixed by the matching conditions. The corresponding term in x 
has the following form 



,2^,2 



5% — c 4 — 2~^ r Ji(V3/cr) cos(2cj t + 25). 
in tt 



H 

This function, as well as its derivatives, is regular for all r including r = 0, and it obeys 
the wave equation 

-$X + $x" $X - k 2 5x = 

r 

in the whole space. Therefore, we shall regard 5x merely as an artifact of the linearity of 
Eq.(57), and interpret 5x as a free, standing electromagnetic wave which is not related to 
the excited vortex. Because we are interested in the excited vortex alone, we put 

c 4 = 0. (74) 
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Thus, we have finally obtained that 



K{r) 



\ 



1 - ^(Mv^r) - I^Ln^kt)). (75) 
m z H 8 



The function h^z) vanishes exponentially for large z, but von Neumann function 
Ni(z) decreases rather slowly for large z, 




mix) = 

Thus, the x field has the long range component equal to 

8 m H m H 

Its natural interpretation is that it describes radiation of the vector field from the excited 
vortex. 

It is easy to check that the resulting total field configuration has infinite energy per 
unit length of the vortex. This is due to the fact that we have considered the excitation 
of the standing wave type for the straight-linear, infinite vortex in infinite space. The 
presence of the long range radiative component means that the excited vortex should not 
be regarded as a localised soliton, eventhough it has the topological charge +1. 



5 Remarks 

Let us recapitulate results of our paper. First, on a methodological side, we have shown 
that with the help of the polynomial approximation one can obtain analytic description 
of the excited vortex, including the profile and frequency of the excitation as well as the 
back-reaction of the excitation on the vortex. Thus, the polynomial approximation turns 
out to be a useful tool also for investigations of dynamics of vortices - earlier it has been 
applied to domain walls 



Second, on a physics side, we have shown that the excited vortex contains the radia- 
tive component. To find it, one has to calculate the back-reaction of the excitation on the 
vortex. The back-reaction is due to non-linearity of Eqs. (26-28). In the linear approxima- 
tion considered in Section 3 the effects which are due to the back-reaction are not taken 
into account. The radiative component we have found is formed of the vector field. We 
expect that after several iterations of our procedure, i.e. correcting the vortex fields, then 
calculating corrections to the excitation field A, next again correcting the vortex fields, 
etc., one would also obtain a radiative component in the Higgs field. Its amplitude will 
be proportional to a higher power of the square of the dimensionless amplitude q 2 o%/m 2 H 
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of the excitation. Taking the back-reaction into account has resulted in drastic alteration 
of physical characteristics of the excited vortex obtained in the linear approximation. 

In the present paper we have considered for simplicity only the standing wave type, 
x 3 -independent solutions in the infinite space, and we have found that in this case the ex- 
citation has infinite energy because of the radiative component. Nevertheless, the obtained 
solution might be of physical interest. For example, in a real system like a superconductor 
of finite volume one could sustain the excitation by an external electromagnetic wave and 
the standing wave regime could be obtained in the whole volume of the superconductor. 
It is clear that obtaining the standing wave in the infinite volume would require infinite 
energy. 

For description of the excitation created dynamically during, e.g., vortex reconnections, 
the standing wave, x 3 -independent solution has to be modified. In such an event the 
excitation is created almost suddenly at certain instant of time with a finite amount of 
energy and it occupies only a small part of the vortex. It turns out that it is not difficult 
to appropriately generalize our calculations, so that one can give approximate, analytic 
description of propagation of such a localised excitation along the vortex, as well as of 
propagation of radiation wave emitted from the travelling excitation Solutions of this 
type have finite energy. 
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